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Abstract. We report on the first detection of the Vela pulsar in the near- infrared with the VLT/ISAAC in the J s 
and H bands. The pulsar magnitudes are J B — 22.71 ± 0.10 and H = 22.04 ± 0.16. We compare our results with 
the available multiwavelength data and show that the dereddened phase-averaged optical spectrum of the pulsar 
can be fitted with a power law F v oc v~ av with a„ = 0.12±0.05, assuming the color excess Eb-v = 0.055±0.005 
based on recent spectral fits of the emission of the Vela pulsar and its supernova remnant in X-rays. The negative 
slope of the pulsar spectrum is different from the positive slope observed over a wide optical range in the young 
Crab pulsar spectrum. The near-infrared part of the Vela spectrum appears to have the same slope as the phase- 
averaged spectrum in the high energy X-ray tail, obtained in the 2 — 10 keV range with the RXTE. Both of 
these spectra can be fitted with a single power law suggesting their common origin. Because the phase-averaged 
RXTE spectrum in this range is dominated by the second X-ray peak of the pulsar light curve, coinciding with 
the second main peak of its optical pulse profile, we suggest that this optical peak can be redder than the first 
one. We also detect two faint extended structures in the l'/5 — 3'.'l vicinity of the pulsar, projected on and aligned 
with the south-east jet and the inner arc of the pulsar wind nebula, detected in X-rays with Chandra. We discuss 
their possible association with the nebula. 
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1. Introduction 

The Vela pulsar was firmly identified in the optical range 
by the detection of optical pulsations with the radio pul- 
sar period IjWallace et al. 19771 IManchester et al. 1978|l . 
The identification has been further confirmed by measur- 
ing the proper motion and parallax of the radio pulsar and 
its optical counterpart ( |Legge 2000| ICaraveo et al. 200 l|l . 
by broad-band photometry revealing peculiar colors of the 
counterpart typical for the optical emission of rotation 
powered pulsars IjLasker 19761 |Mignani &: Uaraveo 2 001 
and refs. therein), and by the high polarization of the op- 
tical emission (Wagner & Seif ert 2000| ). The Vela pulsar is 
an intermediate age, ~ 10 4 yr, isolated neutron star (NS) . 
Its parameters are listed in Table □ 

In comparison with older pulsars detected in the op- 
tical range (see, e.g., |Mignani et al. 2000| ), the Vela pul- 
sar with V= 23™6 is brighter by at least 1 — 2 stellar 
magnitudes. However, available spectral information on its 
optical emission has been limited to broad-band UBVRI 
photometry ( |Mignani fc Caraveo 2001] ), which suggests a 
flat optical spectrum typical for young Crab-like pulsars. 



Send offprint requests to: Yuri Shibanov, 
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* Based on observations collected at the European Southern 
Observatory, Paranal, Chile (ESO Programme 66.D-0568). 



This is in contrast to the middle-aged, ~ 10 5 yr, pul- 
sars PSR B0656+14 and Geminga whose broad-band spec- 
tra are less monotonous, and to the strong excess in the 
near infrared (IR) part of the spectrum of PSR B0656+14 
( |Koptsevich et al. 200 1| ). This may be an evidence of spec- 
tral evolution of the optical emission with pulsar age. The 
young Crab pulsar shows no excess in the IR. Thus, the 
extension of the spectrum of the Vela pulsar towards the 
IR is useful to determine whether the optical properties 
of this pulsar are closer to those of younger or middle- 
aged NSs, and to get additional constraints on the pulsar 
spectral evolution with age. 

The Vela pulsar has also been studied in the 
high energy range, from soft X-rays to 7-rays (e.g., 
IPavlov et al. 2001bl |Harding et al. 2002| and refs. 
therein). Its multiwavelength spectrum is presum- 
ably nonthermal with different slopes in different high 
energy ranges. An exception is the soft X-ray range where 
a strong excess over a power law 1 background is believed 
to be due t o the thermal emiss ion from the surface 
of the NS ( |6gelman et al. 1993| IPavlov et al. 2001 bj) . 
Different slopes imply different nonthermal radiation 
mechanisms at work in the magnetosphere of the pulsar 
(e.g., synchrotron, curvature, inverse Compton scat- 



Hereafter PL, F u oc v 
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Table 1. Parameters of the Vela pulsar (PSR B0833— 45; [Taylor et al. 1993| unless specified otherwise). 
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tering radiation, etc.). These mechanisms are involved 
in different ways in the two competing models for the 
non-thermal emission of pulsars, the "polar cap" model 
(e.g., |Daugher ty fc Harding 19 96|) and the "outer gap" 
model (e.g., |Cheng et al. 1986|IRomani 1996jl . 

It is not yet clear which of these competing models best 
represents the observations, and which radiation mecha- 
nisms are actually responsible for the observed emission in 
each band. In this respect, observations in the IR are im- 
portant to get additional constraints on these mechanisms 
and models. For instance, in the polar cap model IR pho- 
tons, as well as optical ones, can be produced by inverse 
Compton scattering of softer photons by primary and/or 
pair cascade relativistic particles in the magnetosphere. 

A symmetrical Crab-like pulsar wind nebula (PWN) 
with a torus and jet structure has recently been de- 
tected in X-rays with the Chandra observatory around 
the Vela pulsar l|Helfand et al. 20011 IPavlov et al. 2001bl 
IPavlov et al. 2001ajl . Being fainter than the Crab PWN, 
the Vela nebula has not yet been detected in the optical 
range, perhaps because the pulsar field is crowded by back- 
ground stars. Some of the structures of the Crab PWN, 
identified in the optical, appear brighter in the near-IR 
range ( |Sollerman fc Flyckt 2002"| ), showing a PL spectrum 
with a negative slope close to that observed in X-rays. 
There are indications that the X-ray Vela PWN has a PL 
spectrum of a similar slope QGotthelf fc Olbert 2002). In 
this context, deep imaging of the Vela pulsar field in the 
near-IR might be more promising than imaging in the vi- 
sual range for detection of the Vela PWN and for studying 
the mechanisms of the interaction between the relativistic 
pulsar wind and the ambient matter. 

Here we report on the first detection of the Vela pulsar 
in the near-IR J s and H bands, obtained with the VLT. 
The observations, data reduction, astrometry, and pho- 
tometry are described in Sect.|21 The results are discussed 
in Sect. [3] in conjunction with the data available in other 
spectral bands, and summarized in Sect.0] 



2. Observations and data analysis 

2.1. ISAAC observations and data reduction 

The field of the Vela pulsar was observed during three 
nights, December 14 and 15, 2000, and January 5, 2001 



with the Infrared Spectrometer And Array Camera 2 
(ISAAC) attached to the Antu Telescope (UT1) of the 
European Southern Observatory's Very Large Telescope 
(ESO VLT). A log of the observations is given in 
Table In the SW (Short Wavelength) imaging mode, 
the Rockwell Hawaii HgCdTe 1024x1024 array detec- 
tor was used. The pixel size on the sky was 0'.'147 and 
the field-of-view was 2'. 5x2'. 5. The observations were per- 
formed in the J s and H bands in jitter mode, with a 
jitter box size of 20". The J s filter is being established 
as the new standard broadband filter at w 1.24/im by 
most major observatories (Keck, Gemini, Subaru, ESO), 
and is photometrically more accurate than the classical 
J because it is not cut off by atmospheric absorption 
dSimons fc Tokunaga 2002| ILabbe et al. 2003|l . The de- 
tector integration times (DITs) were 45 s and 13 s in the 
J s and H bands, respectively. Each observation was built 
up by a number of DITs per exposure (NDITs), where we 
used 4 NDITs for the J s band and 6 NDITs for the H 
band. The number of exposures (NEXP) in each observa- 
tional block (ObsID) was 13 in J s and varied from 12 to 26 
in the H band. Total exposure times (NDIT x DIT x NEXP 
summed over all ObsIDs) were thus 7020 s and 8268 s in 
the J s and H bands, respectively. 



2 See http://www.eso.org/instruments/isaac/ for details on 
the instrument, filters and observational technique. 



Table 2. Log of VLT/ISAAC observations of the Vela 
pulsar. 

Band Date Time" Exposure Airmass Seeing 6 

UT UT s arcsec 

~T S 14.12.00 7:45 2340 1.072 0.7 

15.12.00 6733 2340 L084 06 
7:18 2340 1.069 0.5 

H 15.12.00 8:03 1014 1.076 0.5 
8:29 1092 1.094 0.5 
8:54 936 1.118 0.7 

05.01.01 6:20 2028 1.073 0.5 
7:05 2028 1.103 0.4 

7:51 1170 1.150 0.4 

a Refers to the first image of the ObsID. 

b Full width at half maximum of the stellar profile. 
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Fig. 1. Large panel: 30" x 30" overview of the Vela 
pulsar field extracted from the 2'.5 x 2'5 frame of the 
VLT/ISAAC image in the J s band. The region bounded 
by the 6"-box is enlarged in the small panels. The same 
region in different bands is shown in Figs. and ^ The 
region bounded by the 10"-box is shown in Fig. [3 Small 
panels show the 6" x 6" field of view centered at the pul- 
sar position in the J s and H bands. Images with the pul- 
sar counterpart subtracted are also shown in the lower 
panels for each band, and the corresponding contour 
plots are given. Reference frame for the contour plots 
shows the image scale only. The pulsar and nearby ob- 
jects ol, o2 and o3 are marked in the H band contour 
plot. The crosses in the J s contour plots show the pulsar 
posi tion as derived from the HST/WFPC2/F555W im- 
age l|Caraveo et al. 2001|) with an uncertainty negligible 
for this image scale. The arrow in the J s contour plot 
shows the pulsar proper motion direction, see Sect 12.21 
for details. 
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The data were reduced with the Eclipse 3 and 
NOAO IRAF software. Dark and flatfield images were 
prepared using standard Eclipse recipes. Then each 
ObsID was considered separately. Image offsets were 
determined using the geomap routine and four field 
stars. The sky background level was determined and 
the images were summed using the routine jitter 
with the parameters Re jectHalf Width = 7, RejectMin 
= 2 , Re j ectMax = 4, which were chosen based on the re- 
sulting image statistics. The parameters of the fully re- 
duced images for each ObsID are summarized in Table [2] 
Finally, offsets between these images were determined and 
the images were combined. Parts of the resulting images 4 
containing the pulsar are shown in Fig. ^ Contour maps 
of these images are also presented in this figure. Isophotes 



3 See http://www.eso.org/projects/aot/eclipse/ 

4 Images are available in FITS format at 



|http://www.ioffe.ru/astro/NSG/obs/vela-ir/| 



of the contour maps correspond to the levels (in counts) 
above the background l n = S + na, where S is the mean 
sky value near the pulsar, a is the sky standard deviation 
per pixel, and n — 1, 2, . . . , 6. 



2.2. Astrometry and morphology of the pulsar field 

For astrometrical referencing of the VLT images we used 
the HST/WFPC2 im age obtained on January 15, 2000 
l|Caraveo et al. 200l)l . The pulsar is clearly detected in 
this image. Positions of 11 reference stars from the im- 
age were used to construct the coordinate transformation 
between the HST and VLT images with the IRAF routines 
geomap/geoxytran. The rms errors of the transformation 
were less than one third of the ISAAC pixel size in both 
RA and Dec. The pulsar position in the ISAAC images at 
the epoch of the VLT observations was calculated using 
the pulsar pixel coordinates in the HST image and the 
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Fig. 2. 10" x 10" vicinity of the Vela pulsar in the J s band. 
Objects o2, o3, and four stars in the upper left corner of 
the image are subtracted (cf. Fig. 1). The arrow shows 
the pulsar proper motion direction. A faint thin elongated 
object oriented approximately perpendicular to the proper 
motion is seen within the region bounded by the ellipse. 
It coincides with the central part of the inner arc of the 
PWN detected in X-rays (see Sect. 13.41 for details). 

pulsar proper motion (Caravc o et al. 2001(1 . This position 
is marked by a cross in the contour plots of the J s im- 
ages in Fig. 2] The pulsar counterpart is clearly detected 
with the offsets -0"01(5) 5 and 0"01(2) in RA and Dec, re- 
spectively, from this position. The errors account for the 
centering uncertainties in the ISAAC images, coordinate 
transformation and the pulsar proper motion uncertain- 
ties. 

Three point-like objects, ol, o2, and o3, are detected in 
the 6"x6" vicinity of the pulsar. They are marked in the 
H-band contour plot in Fig. and their offsets from the 
pulsar position are listed in Table |3| The faintest object 
ol is detected only in the H band. 

An extended structure is seen in the 1"5 vicinity south 
of the pulsar. The structure is more pronounced in the H 
band, but can also be marginally detected in the J s image, 
although with a different shape. The examination of the 
structure in each ObsID shows that its shape varies from 
one ObsID to another. To discriminate the pulsar from 
the extended structure, the IRAF/DA0PH0T software was 
used to construct a point spread function (PSF) in both 
bands using 10 field stars. The pulsar spatial profile was 
fitted with this PSF and subtracted from the images. The 
subtracted images are presented in Fig. ^ 

The extension of the structure in the J s band is 
aligned with the south-east counter-jet from the pul- 
sar detected in X-rays with the Chandra observa- 
tory l|Helfand et al. 20011 IFavlov et al. 2001a|) . The X-ray 
counter-jet is directed opposite to the pulsar proper mo- 
tion marked by an arrow in the J s contour plot in Fig.^ It 
extends up to ~ 100" from the pulsar. The ISAAC struc- 
ture can be a near-IR signature of the X-ray jet in the 2" 

5 Hereafter the numbers in parentheses are uncertainties 
referring to the last significant digit quoted, for example, 
0.01(5) = 0.01 ± 0.05, 22.04(18) = 22.04 ± 0.18. 



vicinity of the pulsar. However, visual inspection of the 
2" vicinity of the pulsar in the Chandra/HRC image does 
not reveal any such structure. It may be hidden in the 
complicated pulsar PSF profile of the HRC image. 

We found in the ISAAC images also a hint of a faint 
thin elongated structure, overlapping with the central part 
of the inner arc of the Vela PWN in the Chandra/ACIS 
images IjPavlov et al. 2001af) . The structure is seen within 
the ellipse in Fig. It is aligned approximately perpen- 
dicular to the pulsar proper motion direction, as is the 
X-ray arc, and its offset from the pulsar is 3"1 along this 
direction. The structure is detected at only ~ 2<r level and 
only in the J s band (see Sect. 12.31 for details). However, 
inspection of each separate ObsID image (see Table EJ 
shows that the structure is absent only in the first J s im- 
age, which has the worst seeing, but it is present in the 
two other images. 

To search for the detected extended structures in the 
images in adjacent bands we examined also the archival 
RI band images of the pulsar field, obtained with the 
VLT/FORS 6 on April 12, 1999 j Wagner fc Seifert 2000| ), 
and in the F675W (overlaps with R) and F814W (over- 
laps with I) bands, obtained with the HST/WFPC2 7 
on March 19, 2000 and on March 15, 2000, respec- 
tively ( Mign anTfc Caraveo 200 1| ) . The reduced images are 
shown in Figs. El and 0] The pulsar is reliably detected 
in all bands, while the extended structures and ol are 
not seen in any of them. The extended structure near 
the pulsar could not be seen even after the pulsar PSF 
was subtracted in the RI bands (Fig. 0J. The object o2 
is seen in both HST bands (integrated exposure time 
2600 s) , but it is only barely visible in the VLT I band and 
not detected in the short R band exposure (300 s). This 
object and the extended structures are also not visible 
in the HST/WFPC2/F555W image llCaraveo et al. 200111 . 
This means that the detected extended structures, as 
well as ol, are red objects. To conclude whether these 
objects are associated with the pulsar nebula or they 
are background objects, additional observations are re- 
quired. A change in brightness of these objects would 
strongly support their association with the highly vari- 
able structure of the PWN, as has been observed in X- 
rays l|Pavlov et al. 2001afl . Note that wisp structures have 
been detected and studied close to the Crab pulsar in X- 
rays and in the optical plester et al. 200*2)1 . and in the IR 
flSollerman fc Flyckt 2002| ). In X -rays and in the optical, 
the observations show that the wisps vary in flux on a 
time scale of about one day. 



6 Based on ESO programme 63.P-0002. 

7 Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained from the data archive at the Space 
Telescope Institute. STScI is operated by the association of 
Universities for Research in Astronomy, Inc. under the NASA 
contract NAS 5-26555. 
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Fig. 3. 6" x 6" vicinity of the pulsar, as in small panels of 
Fig. □ obtained with th e VLT/F0RS1 in R and I bands 
(Wa gner fc Seifert 2000| ). Left panels show original im- 
ages. Right panels show images with the modeled pulsar 
profile subtracted. 



F675W 




Fig. 4. 6" x 6" vicinity of the pulsar, as in small panels of 
Figs.E and in Fig. El obtaine d with the HST/WFPC2 in 
the F675W and F814W bands ( |Mignani fc Caraveo 2001) ). 



2.3. Photometry 

The photometric zeropoints for the ISAAC observations, 
J s0 = 24.81(4) and H Q = 24.56(7), were derived us- 
ing images of the standard star FS13, observed on 
December 15, 2000. The J band catalog magnitude of 
FS13 (Hawa rden et al. 2000j) was used for the J s band. 
Differences in zeropoints between different ObsIDs were 
estimated comparing magnitudes of 6 field stars. The dif- 
ferences proved to be no larger than 1% in the J s and 
2% in the H bands, and were accounted for only in the 
resulting photometric error budget. The average Paranal 
atmospheric extinction of 0.06 mag airmass -1 was used in 
both bands 8 . The fluxes of point-like objects were mea- 
sured in apertures with diameters closest to mean see- 
ing values (see Table El , i.e., 4 pix in J s and 3 pix in H. 



Aperture corrections for the magnitudes were determined 
using 3 field stars. These stars, as well as the aperture 
diameters, were chosen to minimize the uncertainties of 
aperture corrected magnitudes of the faint pulsar, in both 
filters simultaneously. 

An additional check of the photometry was performed 
using the IR-survey 2MASS 9 . Since the region of the 
Vela pulsar is not yet released in the archive catalog, 
we performed photometry of 10 stars present in the J 
and H bands of the 2MASS "Quicklook Images" and in 
our ISAAC images. Their comparison shows that magni- 
tude discrepancies are ( J s )isaac — ^2MASS = 0.15(12) and 
#isaac — #2MASS = 0.03(10). The discrepancy in the J 
band can be partly attributed to the difference in through- 
put of the J and J s filters. Although photometry on the 
2MASS "Quicklook Images" is not recommended because 
of their poor quality, we regard the obtained coincidence 
at the la level as a confirmation of the accuracy of our 
photometric referencing of the ISAAC data. 

The pulsar magnitudes, determined with aperture pho- 
tometry, are J s = 22.61(9) and H = 21.90(13). In addi- 
tion, the pulsar profile was fitted with the PSF constructed 
using IRAF/DA0PH0T. The PSF-fitted magnitudes are J s = 
22.71(10) and H = 22.04(16). They differ by la from 
the aperture photometry magnitudes. This difference may 
be attributed to the contamination of the pulsar counts 
measured with the aperture photometry by the extended 
structure near the pulsar. We therefore consider the PSF- 
fitted magnitudes more reliable. Using these magnitudes 



and the flux zeropoints by van der Bliek et al. (1996) the 
measured fluxes from the pulsar are Fj s = 1.39(12) /j,3y 
and Fh = 1.64(25) /zJy. We performed aperture photome- 
try of the nearby objects ol, o2, and o3, which are marked 
in Fig. ^ The results of the photometry are summarized 
in Table El 

We also measured the surface brightness of the ex- 
tended source in the l'/5 vicinity of the pulsar. It was mea- 
sured on the images with the pulsar subtracted over the 
area (2.4 arcsec 2 , the same in both bands) which covers 
the brightest parts of the structure. The surface bright- 
nesses are 22.98(5) mag arcsec -2 or 1.08(5) /xJy arcsec -2 , 
and 21.51(7) mag arcsec -2 or 2.66(16) /xJy arcsec -2 , in 
the J s and H bands, respectively. The surface brightness 
of the second structure projected at the PWN inner arc 
(see Fig. EJ in the J s band is 25.8(7) mag arcsec -2 or 
0.08(4) /iJy arcsec -2 , and the respective upper limit in 
the H band is 24.44 mag arcsec -2 or 0.18 /_Jy arcsec -2 . 

Since no photometric standards were observed during 
the night of the VLT observations in the RI bands, the 
photometric equations were determined from the Landolt 
standards 10 , observed on the night before: 



R-r = 26.583(09) 
I — i= 25.664(13) 



0.052(17)(r- i) 
0.068(23)(r-£) 



9 
10 



See http://www.eso.org/instruments/isaac/imaging_stan dardihhnhiMt 1992 ). 



http:/ /irsa. ipac.caltech.edu/ 
Fields of PG1323-085, PG1633+099, 



and PG1657+078 
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Table 3. Photometry of the Vela pulsar (PSR) and the nearby objects ol, o2, and o3, marked in Fig. 2] The offsets 
of the objects from the pulsar position are given in the second column. Each spectral band cell for the PSR and 
ol consists of two pairs of magnitude/flux values, divided by lines: upper pair are measured values, lower pair are 
dereddened values. Dereddening was performed with E B _y = 0.055(5). Each pair consists of the magnitude (upper 
value) and the flux in /xJy (lower value). Only measured magnitudes/fluxes are presented for o2 and o3. All magnitudes 
are measured via aperture photometry, except for the pulsar magnitudes in the J s and H bands, which are measured 
with PSF fitting (see Sect. 12.31 for details). No magnitudes for HST bands are calculated. Empty cells mean that the 
object is not measurable in this band. 



Object 


Offset 


H 


Js 


I 


HST/F814W 


R 


HST/F675W 


PSR 




22.04(18) 
1.64(25) 

22.00(18) 
1.69(26) 


22.71(10) 
1.39(12) 

22.66(10) 
1.45(13) 


22.95(13) 
1.57(18) 

22.86(13) 
1.71(19) 


1.367(34) 
1.497(38) 


23.42(11) 
1.29(12) 

23.29(11) 
1.46(14) 


1.279(36) 
1.456(45) 


ol 


0'.'31 W 
1'.'19 N 


22.64(22) 
0.94(17) 

22.61(22) 
0.97(18) 


>24.10 
<0.386 
>24.05 
<0.403 




<0.131 
<0.144 




<0.130 
<0.148 


o2 


l'/18 W 
0'.'29 N 


20.81(11) 
5.07(47) 


21.94(7) 
2.83(18) 




0.685(35) 




0.241(20) 


o3 


1'.'14 W 
1767 S 


18.85(10) 
30.8(2.6) 


19.92(5) 
18.17(87) 


21.53(4) 
5.84(21) 


5.361(59) 


22.70(7) 
2.52(16) 


2.522(48) 



Here R and / are the Cousins magnitudes, r and i are the 
instrumental magnitudes. The ambient conditions moni- 
tor 11 shows that the average atmospheric extinction coef- 
ficient during the observations of the standards was in the 
range 0.140—0.155 mag/airmass, and at the time of the 
pulsar field observations it was 0.150 mag/airmass. Since 
the difference between them is negligible compared to the 
uncertainties of the equations presented above, we used 
these equations for the photometric referencing of our ob- 
servations without corrections for the extinction variations 
between the nights. The aperture corrections were done 
using a PSF constructed from 4 field stars. The measured 
pulsar magnitudes, R = 23.46(11) and / = 22.90(13), cor- 
respond to the fluxes F R = 1.24(12) /iJy and Fi = 1.65(19) 
/ijy, using the magnitude-flux conversion zero-points pro- 
vided by Fukugita et al. (1995) The flux/magnitudes of 
the pulsar and o3 are shown in Table EH 

Pipeline-provided zeropoints and pivot wavelengths 
were used for the flux calibration of the HST observa- 
tions (2.51 x 10~ 18 erg cm" 2 s _1 A" 1 / 7995 A and 
2.90 x 10" 18 erg cmT 2 s^ 1 A" 1 / 6717 A in the F814W 
and F675W bands, respectively). Aperture photometry 
was performed for the pulsar and the objects o2 and o3. 
Aperture corrections were derived from a nearby rela- 
tively bright star. The measured pulsar flux (see Table 0J) 
in the F814W band is consistent with the published one 
UMignani fc Caraveo 200l| |, while the flux in the F675W 
band is apparently ~25% higher. Both fluxes are compat- 
ible with the less accurate VLT fluxes in the RI bands de- 
scribed above. We measured also the fluxes of the nearby 
objects o2 and o3 and estimated 3cr upper limits of the 
object ol in the HST bands. 



11 Available at http://archive.eso.org/asm/ambient-server 



In Table [3] we also present dereddened magnitudes 
and fluxes for the pulsar and ol using Eb-v = 0.055(5) 
(Ay 0.18, R = 3.1). This corresponds to the column 
density Ah = 3.3(3) x 10 20 cm -2 , derived from the com- 
bined PL + NS-atmosphere spectral fit of the Vela pul- 
sar X-ray data obtained with the Chandra observatory 
(Pa vlov et al. 2001b|l . The extinction value is consistent 
with the new distance to the Vela supernova remnant 
(SNR) of 250 ±30 pc l|Cha et al. 1999|> . and with the high- 
est value over the Vela SNR JV H ~ 6 x 10 20 cm" 2 (Ay ~ 
0.32) found for its southern part IjLu & Aschenbach 2000). 

The colors of the brightest stellar object in the pulsar 
vicinity, o3, suggest that it could be a main sequence K5- 
K7 star at a distance of ~ 9 kpc, assuming Ay w 2. This 
is a much higher extinction than we have adopted for the 
pulsar, but it is consistent with the maximum possible 
Galactic extinction in the Vela direction, Ay w 4 (Ah — 
7.5 x 10 21 cm~ 2 ; e.g., |Schlegel et al. 1998| ). The colors of 
the fainter object o2 are roughly consistent with a cooler 
and even more distant main sequence star (of spectral type 
M0) at ~ 10 - 11 kpc (Ay « 3). The object ol is too red 
to be consistent with any ordinary Galactic star. It may 
be associated with the pulsar nebula. It may also be a 
background extragalactic object, as could also the objects 
o2 and o3. We discuss this further in Sect. 13.41 



3. Discussion 

3.1. Multiwavelength spectrum of the Vela pulsar 

In Fig. we have combined our IR data with the available 
phase-averaged multiwavelength fluxes of the Vela pulsar 
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Fig. 5. Dereddened multiwavelength spectrum of the Vela pulsar obtained with different telescopes as marked in the 
plot. Diamond-shaped fillings represent la confidence regions of the NS-atmosphere + PL fit (dot-dashed line) of the 
Chandra data; double-dot-dashed and dotted lines show the contributions of nonthermal PL and thermal atmosphere 
components, respectively. The dashed line shows the PL fit of the near-IR and RXTE data and the solid line represents 
the sum of this PL component with the atmosphere component. The optical and RXTE ranges are shown enlarged in 
Fig. El 



including the radio (EPN 12 ; S. Johnston, 2002, private 
communications), the optical (Mig nani &: Caraveo 20 01), 
X-rays from the Chandra l(plrvToT^t™aT^00TbjL hard 
X-rays from the RXTE fHarding et al. 2002| ) and 
OSSE UStrickman et al. 1996|), and 7-rays from the 
COMPTEL HSchonfelder et al. 2000jl and the EGRET 
( Kanbac h et al. 1994J) . Here we present unabsorbed fluxes 
dereddened with the color excess Eb-v = 0.055(5) ap- 
plied to both the optical and X-ray regions. The VLT 
data in the RI bands are omitted since they are compat- 
ible with the more accurate HST fluxes in the respective 
F675W and F814W bands (cf. TableEJ). 

While the pulsar flux generally decreases with increas- 
ing frequency, one can resolve several nonthermal spec- 
tral components with different slopes in different spec- 
tral ranges. They are presumably of the pulsar magne- 
tospheric origin. An excess in soft X-rays is attributed 
to the thermal emis sion from the surface of the NS 
(Ogelman et al. 1993|lPavlov et al. 2001bjl. 



12 European Pulsar Network (EPN) archive available at 
http:/ /www. mpifr-bonn.mpg.de/div/pulsar/data 



3.2. Phase averaged spectra of the pulsar in the 
optical and X-rays 

As was noted by Pavlov et al. (2001b), Mignani & Caraveo 
(2001), Harding et al. (2002), and seen from Fig. the 
optical emission of the Vela pulsar is likely to be of non- 
thermal origin and the optical data are roughly compati- 
ble with the low energy extension of the X-ray PL spectral 
component dominating in the 2 — 10 keV range. This may 
suggest a similar nature of the optical and the high energy 
X-ray emission. However, as seen from Fig. [5] because of 
the limited statistics of the available Chandra data in the 
2 — 10 keV range, the extension of the PL component in- 
ferred from the Chandra X-ray fit is much less certain than 
the near-IR and optical data. 

At the same time, the RXTE data, being compati- 
ble with the Chandra results, appear to be less uncer- 
tain. A difference in normalization of the PL compo- 
nents detected with Chandra and RXTE can be seen in 
Fig. [5j It may be due to the fact that only the pulsed 
component is detected by RXTE plarding et al. 2002| ). 
But the Chandra observations, representing the total flux 
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Fig. 6. Blow-up of the optical and X-ray parts of the Vela pulsar spectrum presented in Fig. [5J The scale in both 
panels is the same. Left panel: IR-UV part of the pulsar spectrum. IR and optical bands are marked. Diamond- 
shaped, stripe-shaped and solid fillings represent la confidence regions of the NS-atmosphere + PL, IR-RXTE, and 
IR-Optical fits, respectively. The best fits and contributions of different spectral components are shown by different 
types of lines explained in the right panel. Right panel: X-ray part of the spectrum. Fluxes in the RXTE bands are 
marked by crosses. The flux in the hardest RXTE band likely belongs to a flat spectral component detected with the 
OSSE/COMPTEL and was not used in the IR-RXTE fit. Line types are the same as in the left panel and Fig. [3] 



from the pulsar including an "off- pulse" component, show 
that the pulsed fraction in the RXTE range can be as 
high as 80% IjSanwal et al. 200Tj> and the contribution 
of the off-pulsed component may not strongly affect the 
spectral shape. Excluding the 5th most energetic RXTE 
band, which is likely related to a flatter spectral com- 
ponent dominating in the OSSE and COMPTEL ranges 
(see Fig. EJ, the RXTE spectrum can be fitted with a 
PL with the spectral index o£ XTE ~ 0.41(9). Although 
this fit is statistically inconsistent (\ 2 = 3.7 per dof), it 
reflects the spectral slope in the RXTE range: the best 
fit line lies within a narrow stripe-shaped region shown 
at the right panel of Fig. Its extension to the opti- 
cal range overlaps with the optical spectrum. A similar 
behavior was observed for the middle-aged pulsar PSR 
B0656+14 dKoptsevich et al. 200 l\ . But the slope of the 
optical broad-band spectrum as a whole appears to be sig- 
nificantly flatter, a° pt = 0.12(5) (x 2 = 0.7 per dof, see the 
left panel of Fig.EJ, and it is not possible to fit all optical 
and RXTE data with a single PL. We note also deviations 



from the single power law IR-Optical fit at about la level 
seen in the J s and U bands. 

Such behavior of the phase-averaged optical spectrum 
suggests that it can be a combination of several spec- 
tral components dominating at different phases of the 
pulsar light curve, as it is seen in the RXTE range 
(Harding et al. 2002| ). This can only be proven by deep 
time-resolved photometry To our knowledge, no such 
data have been obtained yet for the Vela pulsar. The 
most recent "white-light" time-resolved photometry re- 
veals three peaks in the pulsar light curve in the opti- 
cal range l|Gouiffes 19 98). In contrast to that, up to 5 
peaks were registered in the RXTE bands, and their PL 
spectra have significantly different indices and intensities. 
The second RXTE peak consists of two components, soft 
and hard, and coincides with the second optical peak 
( |Harding et al. 2002| ). The second hard peak dominates 
the whole phase-averaged spectral flux, except for the 5th 
RXTE band, where the first peak with a positive spectral 
slope contributes significantly, providing a smooth connec- 
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tion to the spectral data in the OSSE range (cf. Fig. EJ- 
The first RXTE peak coincides with the first 7-ray peak. 

The measured fluxes in the J s and H bands may imply 
that the pulsar spectrum could be steeper in the IR than in 
the optical, as was also observed for the middle-aged pul- 
sar PSR B0656+I4 ( |Koptsevich et al. 200 1J . Deeper ob- 
servations of Vela in the H band are needed to state this 
possible similarity with greater confidence. What is more 
obvious is that the spectral slope in the near-lR is compat- 
ible with the RXTE slope. Combining the first four RXTE 
bands with the J s and H bands gives a consistent PL fit 
with <4 R " RXTE = 0.417(6) (x 2 = 1.3 per dof) shown by 
a dashed line in Figs. El and El The optical bands show a 
flux excess over this fit. 

To better match the whole optical range we com- 
bined this nonthermal component with the thermal 
NS atmosphere component describing the soft X-ray 
part of the Chandra data (jPavlov et al. 200 lbl combined 
HRC+ASIC fit). The combined model spectrum is shown 
by solid lines in Figs. E] and El However, the atmosphere 
component decreases the residuals significantly only in the 
U and the softest RXTE bands, by contributions from 
the Rayleigh-Jeans and Wien tails of the thermal emis- 
sion, respectively. The rest of the optical bands still show 
a significant excess over the combined atmosphere + PL 
model. Based on that, we can speculate that the 2nd peak 
is responsible for the phase averaged emission in both the 
RXTE and near-lR ranges, while the excess in BVRI is 
mainly produced by another phase and/or spectral com- 
ponent. 

Deeper Chandra observations of the Vela pulsar in the 
high energy tail of its X-ray emission are needed to per- 
form more accurate phase averaged spectral analysis of 
the optical and X-ray data. The high spatial resolution 
of Chandra should avoid the possible uncertainty of the 
RXTE fluxes which do not properly account for the off- 
pulsed component of the pulsar emission. In this context, 
we can assume that the contribution of this component 
just increases the RXTE fluxes by a factor of ~ 1.6 (+0.2 
in Log scale) in all bands to match the Chandra best PL 
fit line (double-dot-dashes in the right panel of Fig. EJ . In 
this case, the low energy extension of the RXTE PL com- 
ponent would only overlap with the upper part of the B 
band error-bar, and would suggest a spectral break of the 
nonthermal component near the UB bands. On the other 
hand, if the current Chandra best PL fit is closer to real- 
ity, the break between the optical and X-ray slopes of the 
nonthermal component would be near 50 eV, i.e., in the 
EUV range. We consider these as alternative hypotheses 
to be tested by future observations. 

3.3. Optical spectrum of the Vela pulsar and spectral 
evolution of the pulsar optical emission 

The near-IR data extend significantly the broad-band op- 
tical spectrum of the Vela pulsar towards longer wave- 
lengths. This allows a detailed comparison with the prop- 
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Fig. 7. Comparison of the optical spectra of four pulsars. 
The youngest Crab ifSollcrma n 2003|l is at the top, the 
oldest Geminga ( Koma rova et al. 2002(1 is at the bottom. 



erties of other pulsars observed in the near-IR. In Fig. [7\ 
we compare the optical-IR part of the spectrum of the 
Vela pulsar with the available phase-averaged optical and 
near-IR spectra of other pulsars of different ages. 

Our photometry of the Vela pulsar in the HST/F675W 
and VLT/R bands does not confirm a dip in this range 
(Mign ani fc Caraveo 2001] ). Moreover, changing the ex- 
tinction from Ay = 0.4 l|Manchester et al. 1978(1 to 
the most recent value Ay = 0.18 changes the aver- 
age spectral index in the optical from ot v = —0.2(2) 
KMignani fc Caraveo 200l\ to a° pt = 0.12(5). This is dif- 
ferent from the most recently estimated positive slope of 
the spectrum of the younger (r w 10 3 yr) Crab pulsar, 
which shows a monotonous flux increase from the IR to 
the FUV range l(Sollerman et al. 20001 TSollerman 2003)) . 
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Spectra of the middle-aged pulsars PSR B0656+14 and 
Geminga have dips near around UB bands, and the spec- 
trum of PSR B0656+14 has a strong increase in the near- 
IR. We see signs of these features also in the Vela spec- 
trum, although at a low significance level. The spectrum 
of PSR B0656+14 increases significantly towards the IR 
(Ko ptsevich et al. 2001| ), contrary to the spectrum of the 
Crab pulsar. 

3.4. Extended structures 

Our IR observations allowed us to detect faint extended 
structures ~l'/5 SE and ~ 3"1 NW of the pulsar. These 
are projected on the SE counter-jet and on the inner 
arc of the PWN de tected in X-rays l|Helfand et al. 20011 
IPavlov et al. 2001a|l . We also see a point-like object ol 
~l'/2 NW of the pulsar, projected on the NW X-ray PWN 
jet. All these objects are red and have no reliable coun- 
terparts in the optical bands. It is difficult to conclude 
whether these features are associated with the PWN or 
they are background objects. 

Observations of the much more energetic Crab PWN 
show that some of its structures are brighter in the IR 
than in the optical range. For instance, the knot structure, 
which is only 0'.'6 SE of the Crab pulsar, is more luminous 
in the IR and has a much steeper spectrum {a^ not ~ 0.8) 
than the Crab pulsar itself ( So llerman 2003f) . This is con- 
sistent with the spatially averaged PL X-ray spectrum 
of the Crab PWN l|Gotthelf fc Olbert 2002jl . Some of the 
Crab wisps are also better resolved in the IR. Based on 
that, and on the faintness of the Vela PWN as compared 
with the Crab one, it would not be too surprising if the 
Vela PWN could be detected more easily in the near-IR 
than in the optical. 

In this context, the measured flux in the H band, 
and our 3a detection limits in the J s , HST/F814W and 
ffST/F675W bands of the point-like object ol resembling 
the Crab knot, suggest a° x ^1.8. Thus, its spectrum is 
much steeper than that of the Crab knot and also steeper 
than a spatially averaged X-ray spectrum of the Vela 
PWN with c^ WN ~ 0.5 HGotthelf fc Olbert 200211 . We are 
not aware of any reasonable physical mechanism for such 
a strong change of the spectral index from X-rays to the 
optical range, and we therefore believe that ol is likely to 
be a distant highly absorbed extragalactic object. 

On the contrary, as seen from Fig.|Hl the surface bright- 
ness in J s and its upper limit in H of the IR inner arc 
shown in Fig.|21are compatible with the low energy exten- 
sion of the time and spatially averaged X-ray spectrum 
of the Vela-PWN inner and outer arcs. The X-ray arc 
spectrum can be fitted with a PL with a u ~ (0.3 — 0.5) 
and unabsorbed energy integrated brightness ~ (3 — 6) x 
10 -14 ergs s _1 cm~ 2 arcscc~ 2 in the 0.1 — 10 keV range 
( |Kargaltsev et al. 2002||Mignani et al. 2003| ). We consider 
this, and the positional coincidence of the inner IR and X- 
ray arcs as arguments in favor of the first detection of a 
counterpart of the Vela PWN in the near-IR range. The 
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Fig. 8. Time and spatially averaged unabsorbed spec- 
trum of the surface brightness of the Vela PWN 
inner and outer arc regions in X-rays fitted by 
a PL ( |KjaTgjdt sev et al. 2002| ) together with optical 
( Mi gnaniTtar^003| ) and radio IjLewis et al. 2002|l upper 
limits, and the brightness of the suggested inner arc coun- 
terpart in the J s band and its upper limit in the H band 
marked by a box. Filled regions show uncertainties of the 
X-ray fit (dashed line) and its extension into the optical 
range. The ellipse outlines the brightness of the near-IR 
extended structure projected at the SE X-ray counter-jet 
in ~l'/5 vicinity of the Vela-pulsar. 



IR brightness of the structure is also consistent with the 
deepest optical upper limit of 0.057 fi3y obtained recently 
in the HST/F555W band ( |Mignani et al. 2003| ). If it is a 
real counterpart of the inner arc with the spectral slope 
described above, only slightly deeper observations, pre- 
sumably at longer wavelengths, would allow a detection 
of this PWN structure. 

The extended IR source closest to the pulsar, and ap- 
parently projected on the SE X-ray counter-jet, is an order 
of magnitude brighter than expected from an extrapola- 
tion of the X-ray spectrum into the near-IR/optical range 
(Fig. |SJ). Its IR spectrum is also much steeper, a v ~ 2.5, 
than the spectrum of the PWN further away from the pul- 
sar. As in the case of ol, this does not argue in favor of 
it being associated with the PWN. It could, however, be 
that the inner jet structure is brighter and has a steeper 
spectrum, because of possible instabilities of the relativis- 
tic particle flow from the pulsar as well as higher radiative 
losses at shorter distances from the pulsar. Although the 
upper limit on the optical flux from the IR source agrees 
with an extrapolation of the flux in the J s and H bands 
(see Fig.|SJl, we emphasize that the optical and near-IR im- 
ages are from different epochs. As mentioned in Sect. 12.21 
the emission of the regions close to the Crab pulsar varies 
on a short time scale, and future comparison between op- 
tical and near-IR emission in the vicinity of the Vela pul- 
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sar would benefit from simultaneous observations in these 
wavelength ranges. 

In X-rays the Vela PWN shows high variability of its 
jet and arc-like structures in position, intensity, and hard- 
ness ratio IjPavlov et al. 2001a)l . Thus, further deep obser- 
vations of the pulsar field in the near-IR and in the optical 
would be useful to search for the variability and to prove 
or reject the association of the detected extended struc- 
tures SE and NW of the pulsar and ol with the PWN. 
Observations of the Vela pulsar in the KL bands would 
be valuable to investigate the possible increase of its flux 
towards the IR range. Finally, time resolved photometry 
and spectral information on the emission of different op- 
tical peaks of the pulsar pulse profile would be crucial to 
understand to which extent nonthermal optical radiation 
of the Vela pulsar is of the same origin as the nonther- 
mal spectral component seen in the high energy tail of its 
X-ray spectrum, or whether it is generated by different 
radiation mechanisms. 

4. Summary 

Here we provide a summary of our most important results. 

1. We have, for the first time, detected the Vela pulsar in 
the near-IR in the J s and H bands. 

2. Our IR fluxes combined with the available broad- 
band optical data confirm the nonthermal origin of 
the pulsar emission in IR-optical range. The combined 
phase-averaged unabsorbed IR-optical spectrum is fit- 
ted with a single PL with a negative slope. This is in 
contrast to a positive slope of the unabsorbed spec- 
trum of the younger Crab pulsar. 

3. The IR-optical spectrum and the phase-averaged PL 
spectral component detected in the high energy tail of 
the pulsar X-ray emission cannot be fitted with a single 
PL. This suggests either a spectral break in the NUV- 
EUV range, or the presence of an additional spectral 
component with a flatter spectrum dominating in the 
optical range. In the latter case, the IR and the X- 
ray spectrum can be fitted with a single PL suggesting 
the same origin of the nonthermal pulsar emission in 
the second pulse of the pulsar pulse profile in both the 
X-ray and IR ranges. 

4. We detected two faint objects in the l'/5 vicinity of the 
pulsar. They are projected on the SE counter-jet and 
the NW jet of the Vela PWN detected in X-rays. Both 
of them are extremely red and have no counterparts 
in the optical range. Their IR fluxes are apparently in- 
consistent with the expected IR brightness of the PWN 
obtained by extrapolation of its X-ray spectrum into 
the IR range. Finding variability of the objects would 
strongly support their association with the highly vari- 
able PWN. 

5. A thin extended structure aligned with the inner arc 
of the X-ray PWN is marginally seen in the J s band. 
Its brightness is consistent with the X-ray PWN spec- 
trum. However, its reality and association with the 



PWN structure need to be confirmed at higher sig- 
nificance level by deeper observations. 
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